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NONCOMMUTATIVE BISPECTRAL DARBOUX 
TRANSFORMATIONS 

JOEL GEIGER, EMIL HOROZOV, AND MILEN YAKIMOV 


Abstract. We prove a general theorem establishing the bispectrality of noncommu- 
tative Darboux transformations. It has a wide range of applications that establish 
bispectrality of such transformations for differential, difference and g-difference op¬ 
erators with values in all noncommutative algebras. All known bispectral Darboux 
transformations are special cases of the theorem. Using the methods of quasideter¬ 
minants and the spectral theory of matrix polynomials, we explicitly classify the set 
of bispectral Darboux transformations from rank one differential operators and Airy 
operators with values in matrix algebras. These sets generalize the classical Calogero- 
Moser spaces and Wilson’s adelic Grassmannian. 


1. Introduction 

1.1. The study of bispectral operators was initiated by Duistermaat and Griinbaum 
in [7] motivated by their applications to computer tomography and time-band limiting. 
Since 1985 the bispectral problem has attracted a lot of attention in pure mathematics 
and has been related to many diverse areas, such as integrable systems (KP hierarchy, 
Sato’s Grassmannian and the Calogero-Moser system) [21 [71 HSl ISH |32], orthogonal 
polynomials [201 |T9], representation theory (VPi+co, Virasoro and Kac-Moody algebras) 
[ana, ideal structure and automophisms of the first Weyl algebra mm and many others. 
Let rii and 02 be two domains in C. The continuous-continuous (scalar) bispectral 
problem asks for finding all analytic functions 

T: Oi X O 2 ^ C 

for which there exist analytic differential operators L{x,dx) and A{z,dz) on Oi and O 2 , 
and two analytic functions 0: Oi —>■ C, / : O 2 — >■ C, such that 

(1.1) L{x,dx)^{x,z) = f{z)^{x,z), 

(1.2) 9{x)'i>{x, z) = A{z,dz)'^{x, z) 

on Oi XO 2 . There are discrete, g-difference and mixed versions of the problem [201 HIES], 
which have played a similarly important role. 

The noncommutative version of the problem is stated for a (generally noncommuta¬ 
tive) associative complex finite dimensional algebra R. It asks for classifying all R-valued 
analytic functions 

T: fli X 02 —^ R 
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for which there exist -R-valued analytic differential operators L{x,dx) and A{z,dz) on 
rii and ^ 2 , and -R-valued analytic functions 0: fli —>■ -R, /: ^2 —>■ -R, such that 

(1.3) L{x, dx)^{x, z) = f{z)^{x, z), 

(1.4) 9{x)'^{x, z) = "^{x, z)A{z,dz) 

on fli X 02- The right action of differential operators in the second equality is defined 

by 

< 6 {x,z) ■ (o(y)S?) := ('['(y. 2)a(2)) 

and linearity. It was observed in [HI EH EH] that the right action was needed in order to 
obtain nontrivial examples which were absent if one insisted on left actions only. More 
precisely, if both actions are left, then they do not commute with each other because 
the algebra R is noncommutative in general. The left action in x and the right action 
in z commute with each other, thus giving rise to a bimodule for two noncommutative 
algebras of differential operators. 

The matrix bispectral problem is the problem in which one takes R = Mn{C) in 
the above definition. It was first considered in [34[ [35] with relation to the ZS-AKNS 
operators. In the last 20 years there has been a great deal of research on the continuous 
and discrete versions of the matrix bispectral problem and numerous relations were found 
with matrix orthogonal polynomials [9ll23[[3n]. spherical functions [22], and many other 
topics. 

The construction and classification of noncommutative (or even just matrix) bispectral 
functions is significantly harder than the scalar problem. One of the major methods 
in the scalar case is the use of Darboux transformations for obtaining large families of 
complicated examples from simpler ones. In this approach, the first equation (II.ip for the 
transformed function is direct while the second one is highly nontrivial. This was carried 
out for second order operators in [7], all rank 1 bispectral operators [31] in Wilson’s 
classification, Bessel functions [2], Airy functions [2[EH]) orthogonal polynomials [20] and 
many other situations. Even though these results are of similar nature, their proofs relied 
on a variety of different algebraic, algebro-geometric and integrable systems techniques. 
A general theorem for bispectrality of scalar Darboux transformations was obtained in 
[T]. In some cases the algebras of dual bispectral operators were also classified [241I261I32] . 

1.2. In this paper we obtain two major sets of results. The first one is a very general 
theorem proving that Darboux transformations from noncommutative bispectral func¬ 
tions produce again noncommutative bispectral functions. All previous examples that 
we are aware of can be recovered as special cases of this theorem, for example all func¬ 
tions in [gilHlEZlISlES]. All results on bispectrality of scalar Darboux transformations 
[H El 0 EQl Ea Ea EH EH also arise as special cases of this theorem. 

Our theorem establishes a general bispectrality type result for all Darboux transfor¬ 
mations in bimodules of noncommutative algebras B. These noncommutative algebras 
B can be taken to be the algebra of analytic differential operators with values in a finite 
dimensional algebra R, the algebras of difference or g-difference operators with values 
in R, and various other examples of a similar nature. For the sake of simplicity in the 
introduction we formulate only the special case of the theorem when the algebra B is 
the first one in the list. We refer the reader to Theorem EH for the general result for 
bispectral Darboux transformations in bimodules of noncommutative algebras. 

Let 'k: Di X D 2 —)■ -R be a noncommutative bispectral function as in (ll.3p - (ll.4p . We 
will assume that it is a nonsplit function in the two variables, in the sense that there 
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are no nonzero analytic differential operators L{x,dx) and A{z,dz) that satisfy ()1.3I) or 
()1.4p with f{z) = 0 or 9{x) = 0. 

To formulate our result, we need to introduce several algebras associated to the func¬ 
tion T(x, z). Let Bi be the algebra of all -R-valued analytic differential operators P{x, dx) 
on rii for which there exits an i?-valued analytic differential operator S{z, dz) on ^2 such 
that 

P{x,dx)^{x,z) = ^{x,z)S{z,dz). 

Let B 2 be the algebra consisting of all such operators S{z,dz)- The nonsplit condition 
implies that the map 

b: Bi B 2 , given by b{P{x,dx)) := S{z,dz), 

is a well defined algebra isomorphism. Denote by M the Bi — B 2 bimodule of i?-valued 
analytic functions on Di x D 2 . Let Ki be the subalgebras of Bi consisting of i?-valued 
analytic functions, i.e. the differential operators in Bi of order 0. Let 

Ai-.= b-\K2), A2:=b{Ki). 

These are the algebras of bispectral operators corresponding to 'I'(x, z). They consist of 
the differential operators that have the properties (ll.3p ~ (ll.4p for some functions f{z), 
9{x). In particular, 

L{x,dx) € Ai and A(z,dz) € A 2 . 

Theorem 1.1. In the above setting, for every factorization, 

L{x,dx) = Q{x,dx)g{x)~^P{x,dx) 

with L{x,dx) G Ai, P{x,dx),Q{x,dx) G Bi, and a regular element g{x) G Ki, the 
R-valued function 

<I>(x, z) := P{x, dx)'It{x, z): Hi X II 2 ^ R 
is bispectral, and, more precisely, satisfies 

{P{x,dx)Q{x,dx)g{x)~^) ^{x,z) = ^{x,z)h{L){z) 

g{x)^{x,z) = ^{x,z) {b{L){z)-^b{Q){z,dz)b{P){z,dz)) ■ 

We recall that a regular element of an algebra is one that is not left or right zero 
divisor. 

The general version of the above result (given in Theorem [H]) deals with arbitrary 
noncommutative algebras Bi,B 2 and Bi — B 2 bimodules M. It presents a general 
approach to establish bispectrality of Darboux transformations in bimodules of non¬ 
commutative algebras. 

Previously, combinations of algebraic, algebro-geometric and integrable systems tech¬ 
niques were used to obtain bispectrality on a case by case basis. This theorem, not 
only proves all these results with a single method, but also has a wide range of new ap¬ 
plications. The theorem directly relates the method of Darboux transformations from 
integrable systems to noncommutative algebra. We expect that it will be also helpful in 
relating bispectrality to matrix factorization |10] . 
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1.3. The second set of results in the paper classifies the bispectral Darboux transfor¬ 
mations in Theorem o in several important cases. The rank 1 bispectral functions of 
Wilson [3T] have played a key role in the classification of scalar bispectral functions. 
They were first introduced in m in terms of wave functions for the KP hierarchy with 
unicursal spectral curves. In [2] it was proved that they are precisely the set of all 
bispectral Darboux transformations from the function Our first result 

classifies all bispectral Darboux transformations from the function 

'k(a;, z) = e^^In in the matrix case R = Mn{C). 

In this case = D 2 = C, and Bi and B 2 are the algebras of matrix differential operators 
with polynomial coefficients in the variables x and 2 , respectively. The isomorphism b 
is given by 

b{d^) = z, b{x) = -d,,b{W) = W, VW€Mn(C). 

The algebras Ki and K 2 are the algebras of matrix-valued polynomials in x and z. The 
algebras Ai and A 2 are the algebras of matrix differential operators in x and z with 
constant coefficients. The classihcation result is as follows: 

Theorem 1.2. The bispectral Darboux transformations from the function e^^In as in 
Theorem \1.1\ for operators L,Q,P with invertible leading terms are precisely the func¬ 
tions of the form 

d>(x,z) = W{x)P{x,dx){e^''In), 

where W{x) is an invertible matrix-valued rational function and P{x,dx) is a monic 
matrix differential operator of order k that has a nondegenerate vector kernel with a 
basis of the form 

for some (not necessarily distinct) aj € C and vector valued polynomials Pj{x) € C[x]"'. 

The operators P{x,dx) are uniquely reconstructed from their vector kernels using 
quasideterminants, see Sect. 13.21 The restriction to operators L, Q, P with invertible 
leading terms is necessary because of the standard pathological problems with leading 
terms that are zero divisors. See Dehnition l3.Il for the notion of nondegenerate space of 
vector-valued meromorphic functions. 

In [311 [32], all rank I scalar bispectral functions were shown to be parametrized by 
the points of the adelic Grassmannian and, equivalently, the points of all Calogero- 
Moser spaces. The classification in Theorem 11.21 and Theorem 4.2 in |33] can be used 
to show that the above class coincides with the recently constructed classes of matrix 
bispectral functions by Wilson m and Bergvelt, Gekhtman and Kasman [5] in relation 
to certain solutions of the multicomponent KP hierarchy and the Gibbons-Hermsen 
system parametrized by the vector adelic Grassmannian. (More precisely, to match 
the two families, our bispectral functions should be multiplied by W{x)~^ and by an 
invertible matrix-valued rational fnnction in z to make the leading terms of their asymp¬ 
totic expansions at 00 equal to e^^In-) One should note that, even though the starting 
function 'I'(x,z) = e^^In is scalar matrix-valued, the set of bispectral matrix Darboux 
transformations from it is much more complicated than the set of bispectral Darboux 
transformation from the latter being the union of all Galogero-Moser spaces |32j . 
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1.4. Our second classification result is for all bispectral matrix Darboux trasnforma- 
tions from the Airy functions. The (generalized) Airy operator is the differential operator 

Af-l 

MAi{x,dx) := ^ aid^~" + aox, 

i=l 

where ao 7 ^ 0,a2, - ■ ■ ,a]sr-i are complex numbers. Its kernel is A^-dimensional over C. 
For any function 'ip{x) in the kernel of MAi{x,dx), define the matrix Airy bispectral 
function 

^Ai{x,z) := 'il){x + z)In. 

The bispectral Darboux transformations from it in the scalar case (n = 1) were classified 
in [21 [28] and played a key role in the treatment of the scalar bispectral problem. Here 
we obtain a classihcation of these transformations for all n. 

In the Airy case, we again have Hi = D 2 = C. The algebras Bi and B 2 are also the 
algebras of matrix-valued differential operators with polynomial coefficients in x and z, 
but the isomorphism 6: Hi —)• i ?2 is given by 

h{x)-.= MAi{z,d,), h{dx)-.= -d,, b{w):=w,yw €MM. 

Here and below MaiIx, dx) is viewed as a scalar matrix-valued differential operator. The 
first two equations imply that b{MAi{x,dx)) = .z. The definition of the isomorphism b 
represents the bispectral equations for the matrix Airy function 

A/Ai(x, (9a;)TAi(x, z) = TAi(x, z)z, xTAi(x, z) = TAi(x, z)Mai{z, dz). 

The algebras Ki and K 2 associated to "^Ai{x-,z) are also the algebras of matrix-valued 
polynomial functions in x and z, respectively. However, in contrast with ^1.31 the 
algebras Ai and A 2 consist of the operators of the form 

q{MAiix,dx)) and q{MAi{z,dz)) 

for matrix-valued polynomials q{t). Denote by 

V’i(x),.. .,'iI^n{x) 

a basis of the kernel of MAi{x,dx)- We have: 

Theorem 1.3. The bispectral Darboux transformations from the matrix Airy function 
TAi(x,z) as in Theorem, \1.1\ for operators L,Q,P with invertible leading terms are ex¬ 
actly the functions of the form 

4>(x, z) = W (x)P(x, 5i,)TAi(x, z), 

where W{x) is an invertible matrix-valued rational function and P{x,dx) is a monic 
matrix differential operator of order dN whose vector kernel is a nondegenerate ndN- 
dimensional space with a basis which is a union of N-tuples of the form 

N-l 

^ A)pj(x), i = l,...N 

j=0 

for some matrix polynomials Pj{x). (The matrix polynomials and complex numbers 
A € C can be different for different N-tuples of basis elements.) 

We expect that the families of bispectral Darboux transformations from the matrix 
Airy functions will have applications to the study of the geometry of quiver varieties 
and to tau-functions of integrable systems satisfying the string equation. 
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1.5. The paper is organized as follows. The proof of Theorem 11.11 and its general 
form for bispectral Darboux transformations in bimodules of noncommutative algebras 
is given in Sect. [2l This section also contains a review of the needed facts from noncom¬ 
mutative algebra. The classification results for matrix rank 1 Darboux transformations 
are in Sect. 01 The classihcation results for bispectral Darboux transformations from 
the matrix Airy function appear in Sect. [5l Sect. [ 6 ] contains examples illustrating the 
classification results. 

The classification results rely on some facts on quasideterminants and the spectral 
theory of matrix polynomials. The former are reviewed in Sect. [3] and the latter in the 
appendix, Sect. [71 

Acknowledgements. We would like to thank Alberto Griinbaum and Alex Kasman for 
their very helful comments on the first version of the paper. M.Y. thanks the hospitality 
of Sofia University, Newcastle University and the Max Planck Institute for Mathematics 
in Bonn where parts of this project were completed. 

2. General Theorem 

In this section we establish a theorem that proves bispectrality of noncommutative 
Darboux transformations in an abstract algebraic setting. After the statement of the 
theorem we explain how it specializes to the concrete settings of continuous and discrete 
matrix-valued bispectral operators. 

2.1. Before we proceed with the statement of the theorem we recall the notion of Ore 
set and noncommutative localization. We refer the reader to m Ch. 6 ] for details. Let 
i? be a unital (in general noncommutative) ring. A subset S of ii is called multiplicative 
if it is closed under multiplication and contains the identity element 1. A multiplicative 
subset E of R consisting of regular elements (non-zero divisors) is called a (left) Ore set 
if for all e £ E and r £ R there exist g £ E and t £ R such that gr = te. If this condition 
is satished, then one can form the localized (or quotient) ring R[E~^] consisting of left 
fractions e“^r, where e £ E, r £ R with a certain natural equivalence relation. The ring 
R canonically embeds in the localization by r i—>■ l~^r. 

For a ring R, we denote by r[R) the multiplicative subset of its regular elements (i.e., 
elements that are not zero divisors). 

Theorem 2.1. Let Bi and B 2 be two associative algebras and M a B\ — B 2 bimodule. 
Assume that there exists an algebra isomorphism b: Bi ^ B 2 and an element £ M 
such that 

PtP = ijb{P), \/P£Bi. 

Assume that Ai and Ki are subalgebras of Bi such that 

b{Ai) = K 2 and b{Ki) = A 2 
and for which the following conditions are satisfied: 

(1) r(Ki) are Ore subsets of Bi consisting of regular elements and 

(2) Ann/(M) = 0 for all f £ r{Ki). 

If L = Qg~^P for some L £ Ai, P,Q £ Bi, and g £ r{Ki), then Lp := Pip = ipb{P) 
satisfies 


{PQg ^)(f = (fb{L), 

gip = <f {biL)-^biQ)b{P)) . 
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Note that b(L) € K 2 . The two conditions ensure that the algebras Bi embed in 
Bi[r{Ki)~^] and M in M[r{Ki)~^][r{K 2 )~^], see e.g. [T71 Ch. 10] for details. Both 
equalities take place in the localized module. 

We prove the theorem in Sect. 12.41 In the next two subsections we describe how it is 
applied in the continuous-continuous and discrete-discrete noncommutative bispectral 
situations. 


2.2. In the continuous-continuous noncommutative bispectral setting the theorem works 
as follows. Let i? be a complex finite dimensional algebra. For i = 1,2, let be two 
open subsets of C. Let 

T: Hi X 02 —-R 

be an i?-valued bispectral function. Recall from the introduction that this means that 
there exist R-valued analytic differential operators L{x,dx) and A{z,dz) on Hi and 0,2, 
and two R-valued analytic functions 

e-.Oi^R, f:Q2^R, 


such that 


(2.1) L{x, dx)'i>{x, z) = T(x, z)f{z) 

(2.2) 9{x)'i>{x, z) ='^{x, z)A{z,dz) 

on Oi X 02- We will furthermore assume that ^{x,z) is a nonsplit function of x and z 
in the sense that it satisfies the condition 

(*) there are no nonzero analytic differential operators L{x,dx) and A{z,dz) that 
satisfy one of the above conditions with f{z) = 0 or 6{x) = 0 . 

In this setting, denote by Bi the algebra of all R-valued analytic differential operators 
P{x,dx) on Oi for which there exits an R-valued analytic differential operator S{z,dz) 
on O 2 such that 

P{x,dx)'^{x,z) = '^{x,z)S{z,dz). 

Denote by B 2 the algebra of all such operators S{z,dz)- It follows from the assumption 
(*) that the map b: Bi —>■ B 2 , given by b{P{x,dx)) ■= S{z,dz), is a well defined algebra 
isomorphism. Let M be the Bi — B 2 bimodule consisting of R-valued analytic functions 
on X D 2 . 

Let Ki be the subalgebras of Bi consisting of R-valued analytic functions. Set 

Ai:=b-\K2), A2:=b(Ki). 

These are the algebras of bispectral operators corresponding to 'I'(x, z). They consist of 
the differential operators that have the properties (I2.1ji - (l2.2jl for some functions f{z), 
9{x). In particular, 

L{x,dx) ^ Ai and A{z,dz) ^ A 2 . 

The sets of regular elements r{Ki) consist of those R-valued analytic functions 9{x) and 
f{z) that are nondegenerate (in the sense that their determinants are not identically 
zero on any connected component of Dj). This immediately implies that the condition 
(2) in Theorem 12.II is satisfied in this setting. By standard commutation arguments with 
differential operators one shows that the subsets r(iLj) are Ore subsets of B^. 

Theorem 12.11 implies the following: 


Proposition 2.2. In the above setting, for every factorization, 

L{x,dx) = Q{x,dx)g{x)~^P{x,dx) 
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for L{x,dx) G Ai, P{x,dx),Q{x,dx) € Bi, g{x) G r{Ki), the R-valued function 

z) := P{x, dx)'it{x, z): Qi X il,2 ^ R 
is bispectral, and, more precisely, satisfies 

{P{x,dx)Q{x,dx)g{x)~^) ^{x,z) = ^{x, z)b{L){z) 

g{x)^{x,z) = ^{x,z) {b{L){z)~^b{Q){z, d„)b{P){z, d^)) . 

2.3. In the discrete-discrete noncommutative bispectral setting Theorem 12.II works as 
follows. As before, let i? be a complex finite dimensional algebra. For i = 1, 2, let flj be 
two subsets of C which are invariant under the translation operator 

Ta;: X i-A X + 1 


and its inverse. Denote by M the space of i?-valued functions on Di x Q 2 - An valued 
difference operator on is a finite sum of the form 

J2ckix)T^, 

fcgZ 

where Ck are i?-valued functions. The space M is naturally equipped with the struc¬ 
ture of a bimodule over the algebras of ii-valued difference operators on Di and D2, 
respectively, by setting 

I ^Cfc(x)T^ I • T(x, 2 ) := ^Cfc(x)T(x + k,z), 

^ix,z)- (j2<^k{z)TA :='^Ckiz - k)^{x,z - k). 

VfcgZ / fcgZ 


An i?-valued discrete-discrete bispectral function is by definition an element of M (i.e., 
an i?-valued function) 

T : X D2 —^ i? 

for which there exist i?-valued difference operators L{x,Tx) and A{z,Tz) on Di and Q 2 , 
and i?-valued functions 


such that 


e-.ni^R, f:n2^R, 


(2.3) L(x, Tx)'it{x, z) = T(x, z)f{z) 

(2.4) 0(x)T(x,z) = T(x,z)A(z,r,) 

on Qi X D2- Similarly to the continuous case, we will assume that 'I'(x,2;) is a nonsplit 
function of x and z in the sense that it satishes the condition 

(**) there are no nonzero difference operators L{x,Tx) and A{z,Tz) that satisfy one 
of the above conditions with f{z) = 0 or 9{x) = 0. 

Let Bi be the algebra of all i?-valued difference operators P{x,Tx) on Di for which 
there exits an valued difference operator S{z,Tz) on D2 such that 

P{x,Tx)'i>{x,z) = 'i>{x,z)S{z,Tz). 

Let B 2 be the algebra of all such operators S{z,dz). The assumption (**) implies 
that the map b: Bi —)■ B 2 , given by b{P{x,dx)) ■= S{z,dz) is a well defined algebra 
isomorphism. Denote by Ki be the subalgebras of Bi consisting of i?-valued functions. 
The algebras 


A,-.= b-HK2), A2-.= b{Ki) 
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consist of the bispectral difference operators corresponding to '^{x,z) (i.e., difference 
operators in x and z that have the properties (I2.3p ~ (l2.4p i. The sets r{Ki) of regular 
elements of the algebras Ki consist of the nonvanishing functions 9{x) and f{z). It is 
easy to verify that the second condition in Theorem 12.II is satisfied in this situation and 
that r{Ki) are Ore subsets of Bi. We leave the details to the reader. 

Theorem 12.11 implies the following: 

Proposition 2.3. In the above setting, for every factorization 

L{x,T,;) = Q{x,T,,)g{x)~^P{x,Ta;) 

with L{x,Tx) G Ai, P{x,Tx),Q{x,Tx) G Bi, g{x) G r{Ki), the function 

^{x,z) := P{x,T,c)^{x,z) 

is an R-valued discrete-discrete bispectral function on fli xn 2 - More precisely, it satisfies 

{P{x,T,,)Q{x,T,,)g{x)~^) ^{x,z) = ^{x, z)h{L){z), 

g{x)^{x,z) = ^{x,z) {b{L){z)~^b{Q){z,Tz)b{P){z,Tz)) . 

Analogously, Theorem 12.11 applies to noncommutative bispectral functions for q- 
difference operators and to every mixed situation, where the operators on the two sides 
are differential, difference, or g-difference of different type. We leave to the reader the 
formulations of those results, which are analogous to the ones above. 

The bispectral functions in [6l[l8l|271|3l] are special cases of the continuous-continuous 
version of the theorem and the ones in p n Eu Ea Ea SO] are special cases of the 
discrete-continuous version. 

2.4. We finish this section with the proof of Theorem 12.11 

Proof of Theorem \2.1\ The first statement is straightforward, 

{PQg~^) if = {PQg~^) Pif = PLif = Pfjb{L) = ipb{L). 

The key point of the theorem is the second statement. To prove it, first note that 
since b: Bi ^ B 2 is an isomorphism, the subalgebras Ai of Bi also satisfy the conditions 
(l)-(2). The Ore subset assumptions imply that 

g = PL~^Q and thus b{g) = b{P)b{L)~^b{Q). 

Therefore, 

g^ = gf;b{P) = 'ifb{P)b{L)-H{Q)b{P) = ip {b{L)-^b{Q)h{P)) . 

□ 

3. Quasideterminants, kernels and factorizations of matrix differential 

OPERATORS 

In this section we gather some facts relating the factorizations of matrix-valued dif¬ 
ferential operators and their kernels, which will be needed in the following sections. 
These results are minor variations of those obtained by Etingof, Gelfand and Retakh 
in HU. One difference in our treatment is that we focus on vector kernels rather than 
matrix kernels (considered in [TU) because our classihcation results for bispectral Dar- 
boux transformations from the next 2 sections are naturally formulated in terms of the 
former. 
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3.1. We first recall the notion of quasideterminants introduced and studied by Gelfand 
and Retakh [laiii]. Let R be an associative algebra. Let X = (x^) be a matrix with 
elements in R. For each pair of indices i,j denote by Ti{X) the i-th. row of X and 
by Cj{X) the j-th column of X. Denote by X^^ the matrix obtained by removing the 
f-th row and the j-th column of X. For each row r and index j, dehne the vector 
obtained by removing the j-th entry. Similarly, for a column c, denote by the column 
vector obtained by its i-th entry. Finally, we assume that the matrix is invertible. 
The quasideterminant \X\ij is defined by 

\X\,, := G R. 

Let D C C be a domain and R be the algebra of meromorphic functions on D with 
values in M„(C). 

For any sets of matrix-valued meromorphic functions ,..., on D and vector¬ 


valued meromorphic functions fi,, f^k on D, define the Wronski matrices 
(3.1) 




( F, . 

. Fk \ 



/ fl ■ 

fnk \ 

W{Fu. 

■,Fk) = 

Fi . 

• f;^ 


■ ■ ) fnk) — 

fi ■ 

f 

J nk 




• Ft^^/ 




Ank-l) 

■ Jnk / 


Such a set of functions Fi,... ,Fk is called nondegerate if W{Fi,... ,Fk) is an invertible 
matrix-valued meromorphic function on D. The invertibility condition is independent 
of whether (j3.1h is considered as an n x n matrix with values in R or an nk x nk matrix 
whose entries are meromorphic functions on D. This is equivalent to saying that the 
determinant of the Wronski nk x nk matrix is a nonzero meromorphic function on D. 

Definition 3.1. A subspace of dimension nk of the space of vector-valued meromorphic 
functions on D will be called nondegenerate if it has a basis fi,..., f^k such that the 
Wronski matrix W{fi,..., fnk) is an invertible matrix-valued meromorphic function on 

n. 


It is clear that the condition in the dehnition does not depend on the choice of basis. 

Remark 3.2. For every basis of a nondegenerate subspace of the space of vector-valued 
meromorphic functions on D, there exists an arrangement of the basis elements into the 
columns of a set of meromorphic functions Fi,..., Fk'. D —M„(C) with the property 
that for all m < A: the set Fi,..., Fm is nondegenerate on D. This follows from the fact 
that, for every nondegenerate square matrix M, there exists a permutation matrix S 
such that all leading principal minors of MS are nondegenerate (the fact is applied for 
M equal to the Wronski matrix of the basis elements). 

3.2. Dehne the vector kernel of a differential operator P{x,dx) with matrix-valued 
meromorphic coefficients on D by 

(3.2) vkerP(x,9a;) := {meromorphic functions f{x): D —>■ C"" | P{x,dx)f{x) = 0}. 

By the standard theorem for existence and uniqueness of solutions of ordinary differ¬ 
ential equations, dimvker P(x, Sa,) < kn, and each a € D at which P(x,dx) is regu¬ 
lar has a neighborhood Oa ^ D such that the restriction of P{x, dx) to Oa satishes 
dimvkerP(x,9a;) = kn. 

Proposition 3.3. Let Q CC be a domain and n, k he positive integers. 
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(a) Assume that P{x,dx) is an Mn{C)-valued meromorphic differential operator on 
ri of order k that has nk-dimensional vector kernel. Then its kernel is nondegenerate in 
the sense of Definition \3.1[ 

(b) Assume that V is a nondegenerate subspace of the space of veetor-valued mero¬ 

morphic functions on ri of dimension nk. Then there is a unique monic differential 
operator P{x,dx) of order k with matrix-valued meromorphie eoeffeients on ri whose 
vector kernel equals V . For every basis ofV choose an arrangement of the basis elements 
onto the columns of a set of meromorphic functions Fi,..., Fk- ^ M„(C) with the 

property that for all m < k the set Fi,..., Fm is nondegenerate on Q (cf. Remark \3.Sfl . 
The operator P{x,dx) is given by 

(3.3) P{x,dx)F{x) = \W{Fi,... ,Fk,F)\k+i,k+i 
for each matrix-valued meromorphic function F{x) and by 

(3.4) P{x, dx) = {dx -bk)... {dx - 6i), 
where 

(3.5) bj ■.= W'^W-\ Wj:=\WiFi,..., Fj ) . 

Proof, (a) Denote by the open dense subset of D on which P(x,dx) is analytic. 
Choose a basis /i(x),..., fnkix) of vker P{x, dx). Let cq, ..., Cnk € C. The statement of 
this part of the proposition follows from the fact that if 

W{fi, . . .,fnk){xo) ■ (ci, . . . Cnk)*' = 0 

for some xq € D°, then this holds for all xq G D°, where {.f denotes the transpose of a 
matrix. (This is nothing but the uniqueness statement for solutions of the initial value 
problem for ordinary differential equations.) 

(b) By [11) Theorem 1.1] there exists a unique monic differential operator P{x, dx) of 
order k such that 

P{x,dx)Fi{x) = 0, y 1 < i < k 

and this operator is given by (|3.3p and ()3.4p . Since vkerP DV,nk> dim vker P and 
dim!/ = nk, we have vkerP = V. □ 

3.3. The factorizations of matrix-valued differential operators are described by the next 
proposition. 

Proposition 3.4. Let L{x,dx) be an Mn{C)-valued monic meromorphic differential 
operator on a domain D C C o/ order I whose vector kernel has dimension nl. All 
factorizations 

(3.6) L{x, dx) = Q{x, dx)P{x, dx) 

into Mn{C)-valued monic meromorphic differential operators of orders I — k and k are 
in bijection with the nondegenerate nk-dimensional subspaces of vkerL. This bijection 
is given by P vkerP. The inverse bijection from nondegenerate subspaces o/vkerP 
to monic differential operators P satisfying ()3.6p is given by Proposition \3.A (b). 

This proposition is similar to [161 Theorem 1], but the latter omits the condition for 
nondegeneracy of vker P. 

Proof. Assume that we have a factorization as in (|3.6p . Obviously, vkerP C vkerP. We 
will show that dim vker P = nk. Then by Proposition 13.31 (a), vkerP is nondegenerate. 
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Fix a point a € O at which L{x, dx) is regular and a neighborhood Oa of it on 
which Q and P are analytic and the restriction of P has refc-dimensional kernel. For this 
restriction vker P C vker L, i.e., there exists an n/c-dimensional subspace V of the vector 
kernel of L{x,dx) such that the vector kernel of the restriction of P to Oa is V\oa- By 
applying Proposition 13.31 (b), we construct a differential operator P'{x,dx) on Q with 
meromorphic coefficients such that vkerP' = V. The restrictions of P and P' to Oa 
coincide by Proposition 13.31 (b) and the two differential operators are meromorphic on 
n. Since 17 is connected, P = P'. So, vkerP = V , and, in particular, dim vker P = nk. 

In the other direction, for an nfc-dimensional subspace V of the vector kernel of L on 
fi and a basis of it, we can form meromorphic functions Fi,..., Fk- M„(C) from 
a basis of V as in Remark 13.21 Furthermore, one can complement the basis of P to a 
basis of vkerL and form a set of meromorphic functions Pfc+ij ■ ■ ■ ,Pi- 17 Mn{C) from 
them, such that Pi,... Pm is a nondegenerate set for all m < /. Then, by (13.4p . 

L{x,dx) = Q{x,dx)Pix,dx) 
for 

Q{x, dx) ■■= {dx -k)... {dx - bk+i), P{x, dx) := {dx - bk) ■ ■ ■ {dx - 6i), 
where bj are given by (13.5|) . □ 

In the setting of the above proposition the meromorphic operator PQ is called a 
matrix Darboux transformation from L = QP. For an eigenfuction ^{x,z) of L{x,dx) 
(satisfying L{x,dx)'^{x, z) = '^{x, z)f{z)) the function P{x,dx)^{x, z) is also called a 
Darboux transformation from ^{x,z). 

Remark 3.5. Some authors consider the following form of Darboux transformations. 
For a given matrix-valued monic meromorphic differential operator L{x,dx), consider 
all matrix-valued meromorphic operators L{x,dx) that have the property 

(3.7) P(x, dx)L{x, dx) = L{x, d)P{x, dx) 

for some matrix-valued monic meromorphic operators P{x,dx)- As before, we call the 
function P{x,dx)^{x,z) a Darboux transformation from ^{x,z). 

A priori it appears that this form of Darboux transformations is more general than 
the ones we consider. However, this is not the case, i.e.. Theorem [Q also covers these 
transformations. To see this, one first proves that (I3.7jl implies P(vkerP) C vkerP. 
Since dim vker P < oo, vkerP C vker/i(L) for some polynomial h{z). Therefore h{L) = 
QP for some matrix-valued meromorphic differential operator Q on 17 and we have the 
Darboux transformation 

h{L) = QP ^ PQ = h{L) 

which falls within the class treated in Theorem EH On the level of wave functions it 
represents the Darboux transformation (13.7p 

T(x, z) i-7> P{x, dx)'i’{x, z). 

4. Classification of matrix bispectral Darboux transformations in the 

RANK 1 CASE 

All bispectral algebras of rank 1 scalar ordinary differential operators were classified by 
Wilson m- In [2] it was proved that the corresponding wave functions are bispectral 
Darboux transformations from the function ^{x,z) = in the sense of Theorem 
EH At the time the result was phrased in different terms, but it is easy to convert 
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the statements of [2] to the ones in Theorem 12.11 and more precisely, to the case of 
Proposition 12.21 when the size of the matrices is 1 x 1. 

In this section we classify explicitly all bispectral Darboux transformations of the 
corresponding matrix analog in full generality, i.e., we obtain a full classihcation of all 
matrix rank 1 bispectral Darboux transformations. This is achieved in Theorem 14.31 
The only restriction we impose is that the leading terms of the differential operators L, Q 
and P in Proposition 12.21 be invertible matrix-valued meromorphic functions. Without 
loss of generality this situation is the same as the one when L, Q and P are monic. (The 
classifications of the factorizations of differential operators with noninvertible leading 
terms are related to well known pathological problems from the zero divisors of the 
algebra Mn{C) and will be considered elsewhere.) 

4.1. In the setting of the previous section, let = D 2 := C. Dehne the function 

T(x,2;) = 

on C X C, where In is the identity matrix. Let Bi and B2 be both equal to the algebra 
of matrix-valued differential operators with polynomial coefficients. We will use the 
variables x and 2 for the first and second algebra, respectively. Let b: Bi ^ B 2 he the 
isomorphism given by 

b{d^):=z, b{x) = -d,, b{W):=W,VW €Mn{C). 

We are in the bispectral situation of Proposition 12.21 

P{x,d^)'I){x,z) = ^{x,z){bP){z,dz), yP{x,d^) € Bi. 

Let Ki and K 2 be the algebras of matrix-valued polynomial functions and x and 
respectively. Denote by Ai and A 2 the algebras of matrix-valued differential operators 
(in x and z) with constant coefficients. Then 

b{Ai) = K2, b{K2) = Ai. 

The regular elements of the algebras Ki and K 2 are the matrix-valued polynomial 
functions 9{x) and f{z) that are nondegenerate in the sense that their determinants 
are nonzero polynomials. By standard commutation arguments for matrix differential 
operators the sets r{Ki) are Ore subsets of B^. The localizations Bi[r{Ki)~^] and 
B 2 [r{K 2 )~^] are nothing but the algebras of differential operators with rational coeffi¬ 
cients in the variables x and Theorem 12.11 implies the following: 

Theorem 4.1. If L{x,dx) is a matrix-valued differential operator with constant coeffi¬ 
cients, and Q{x,dx) and P{x,dx) are matrix-valued differential operators with rational 
coefficients, such that 

(4.1) L(x, dx) = Q{x, dx)P{x, dx), 
then 

^{x,z) := P{x,dx){e^''In) 

is a matrix bispectral function. 

More precisely, if 

(4.2) L{x, dx) = Q'{x, dx)g{x)~^P'{x, dx) 

for matrix-valued differential operators with polynomial coefficients Q'{x,dx), P'{x,dx) 
and a nondegenerate matrix-valued polynomial function g{x), then <h(x, z) := P'{x, dx){e^^In) 
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satisfies 

{P'{x,dx)Q'{x,dx)g{x)~^) ^{x,z) = ^{x, z){bL){z), 

g{x)<^>{x,z) = ^{x,z) {{bL){z)~^{hQ'){z,dfi){hP'){z,dfi}) . 

4.2. Our goal below is to classify all wave functions z) that can be obtained from 
Theorem [Q for operators L{x,dx) with nondegenerate leading term. This is done in 
Theorem 14.81 Denote the vector space of vector-valued quasipolynomial functions 

(4.3) QVn ■■= ©„gce“"C[x]’", 

where C[x] denotes the space of polynomials and ©[x]” the space of vector-valued poly¬ 
nomial functions. The elements of QVn have the form e^^pa,i{x), where the sum is 
finite and pa,i{x) G ©[x]”. 

By converting n-th order vector differential equations to first order ones and using 
Proposition 13.41 one easily obtains the following lemma. We leave its proof to the reader. 

Lemma 4.2. (a) Assume that L{x,dx) is a monie matrix-valued differential operator 
with constant coefficients. IfQ{x,dx) o,nd P{x,dx) are monie matrix-valued differential 
operators with meromorphic coefficients in C which satisfy dm), then vkei P{x,dx) is 
a nondegenare subspace of QVn ■ 

(b) For every nondegenerate finite dimensional subspace V of QVn whose dimension 
is a multiple of n, there exist monie matrix-valued differential operators L, Q and P 
satisfying the properties in part (a) such that vker P{x,dx) = V- 

The next theorem is our classification result of the matrix rank 1 bispectral Darboux 
transformations. 

Theorem 4.3. Assume that L{x, dx) is a monie matrix-valued differential operator with 
constant coefficients, and Q{x,dx) and P{x,dx) are monie matrix-valued differential op¬ 
erators whose coefficients are rational functions. If (14.ip is satisfied, then vker P{x,dx) 
is a nondegenerate subspace of QVn that has a basis of the form 

(4.4) {e“'’'*Pj(x) I 1 < j < kn}, 

where ai,..., akn cire not necessarily distinct complex numbers and Pj{x) G ©[x]”. Fur¬ 
thermore, for every such subspace of QVn there exists a triple of operators L, Q and P 
with the above properties such that vkei P{x,dx) equals this subspace. 

Corollary 4.4. The bispectral Darboux transformations from the function as in 

Theorem \4-l\ for operators L,Q,P with nondegenerate leading terms are precisely the 
functions of the form 

4>(x,z) = iy(x)P(x,5a;)(e*^/n), 

where P{x,dx) is a monie matrix differential operator of order k with vector kernel 
having a basis of the form dm and W{x) is a matrix-valued nondegenerate rational 
function. The operator P{x, dx) is given by (13.3p . where Fffx ),..., Tfc(x) are the matrix¬ 
valued functions with columns 

e"i"^pi(x),...,e"'=->fc„(x), 

Oj G C , Pj(x) G CPx]”. 

The bispectral functions in this family were related to wave functions of solutions of 
the n-KP hierarchy by Wilson m and Bergvelt, Gekhtman, and Kasman m- The fact 
that the family in Corollary 14.41 is the same as the one in [5l [33] (up to a multiplication 
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on the left by W(x) ^ and by an invertible matrix-valued rational function in z) follows 
from Theorem 4.2 in m and Corollary 14.41 Since the proof of Theorem 4.2 in [33] is 
only sketched there, we note that the reference to it is not used in proofs here, but only 
to relate the bispectral functions in this section to those in The bispectrality of 

the special elements of this family obtained by transformations from L{x,dx) = q{dx)D 
for a diagonal matrix D and a (scalar) polynomial q{t) was previously obtained in |6|. 
Applications of Theorem 12.11 to generalized bispectrality in relation to the results of 
Kasman m are described in 1)4.41 

4.3. We will break the proof of the Theorem 14.31 into two lemmas. 

Lemma 4.5. Let, for 1 < j < kn, aj € C and pj{x) be vector-valued rational functions 
for which 

(4.5) e“i^pi(x),..., e“'="^pfc„(x) 

form a nondegenerate (and thus, linearly independent) set. Then the coefficients of 
the unique monic matrix-valued differential operator with kernel spanned by (|13|) are 
rational functions. 

Proof. We argue by induction on k. In the case k = 1 the differential operator equals 

dx-F'{x)F-\x), 

where F{x) is the matrix function with columns e“i*pi(x),... ,e°‘'^^pn{x), and (e.g. by 
Cramer’s rule) F'{x)F~^[x) has rational coefficients. 

In the general case, by (|3.4p . the differential operator has the form 

P{x,dx) = P{x,dx){dx - F'{x)F~^{x)) 

with T(x) being nondegenerate and with columns as in ()4.5I) . Therefore, dx—F'{x)F~^{x) 
has rational coefficients and the operator P{x,dx) has vector kernel spanned by 

{dx - F'{x)F~^{x)){e°'i^pj{x)), n < j < kn. 

By the inductive assumption, applied to P{x,dx), the operator P{x,dx) has rational 
coefficients, which implies the same property for P{x,dx). □ 

Lemma 4.6. Assume that P is matrix-valued differential operator with coefficients that 
are rational functions. If 

r 

^^e^^^Psix) € vkerP 

for some vector-valued polynomial functions ps{x) and distinct complex numbers Ai,..., 
then 

e^‘^Ps{x) G vkerP, VI < s < r. 

Proof. We have 

r 

s=i 

Since the functions e~^‘^P{x, dx){e^‘‘^Ps{x)) are rational, the nonzero functions in the 
sum are linearly independent because Ai, A 2 , • • •, A,, are distinct. Thus all terms in the 
sum vanish. □ 


Proof of Theorem 4-3. The first part of the theorem follows from Lemmas 
(a). The second part of the theorem follows from Lemmas 14.51 and 14.21 fbl. 


and 


□ 
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4.4. There are generalized versions of the differential bispectral problem from (II.3p - 
(II.4p . In one of them, one keeps the differential eigenvalue problem in (|1.3p but replaces 
the operator A(z, dz) in the right hand side of (11.41) with a mixed differential-translation 
operator. In this subsection, we briefly describe an application of Theorem 12.II that gives 
a second proof of the results in m on the construction of such generalized bispectral 
functions. 

Firstly, for an analytic function T: Hi x 142 ^ Afjj(C) with an expansion 

OO 

T(x,2)= '^{x-xoyaj^k{z-^ Mn{C),Xo eQi,zo 
j,k=0 

and a matrix W G dehne the translation action 

CXD 

4'(x,z) • (T]^) := ^ (x - xoYttj^kiz - IT - zq)’^ 

j,k=0 

whenever the series converges. A mixed differential-translation operator is an operator 
obtained by composing differential operators with translation operators of this kind. 

For an invertible matrix H G M„(C), denote by Z{H) = {IT G M„(C) \ WH = HW} 
its centralizer. Consider the analytic function 

^h{x,z) 

Let Bi be the algebra of differential operators in x with coefficients that are quasipoly¬ 
nomial functions with values in Z{H) (i.e., functions with values in Z{H) whose matrix 
entries are finite sums of terms of the form x^e°‘^, j G N, a G C). Let B 2 be the algebra 
consisting of mixed differential-translation operators of the form 

a.MdiTf 

jGN,cGC 

given by finite sums and terms aj^dz) that are polynomials with values in Z[H). The 
following define an isomorphism b: Bi ^ B 2 

b{x) = -dz, b{dd = z, 6(e“^) = 6(IT) = IF, VIF G Z{H) 

and the function 'ifnixjz) satishes 

(4.6) P ■'^h{x,z) = ^H -HP), VPgBi. 

Let Ai be the subalgebra of Bi consisting of constant differential operators in x with 
values in Z{H) and K 2 be the subalgebra of B 2 consisting of polynomial in z with values 
in Z{H). The identities in (14.61) for P G Ai are the spectral differential equations in x 

L{x, dd'^nix, z) = z)b{L){z), L(x, d^) G Ai. 

Let A 2 be the subalgebra of B 2 spanned by the differential-translation operators of the 
form diT^^ ^IT, j G N, c G C, IT G Z{H). Let Ki be the subalgebra of Bi consisting 
of quasipolynomials in z with values in Z{H). The identities in ()4.6p for P G Ki are 
the spectral differential-translation equations in z 

6"^(A)(x)Th(x, z) = Th(x, z) • a, a G A 2 . 

We are in the situation of Theorem O and can apply it to obtain functions that 
satisfy spectral differential equations in x and spectral differential-translation equations 
in z by Darboux transformations from 'I'//(x, z). This gives a second proof of the results 
of Kasman m- 
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Theorem 4.7. Let H G M„(C) be an invertible matrix. If L{x,dx) is a Z{H)-valued 
differential operator with constant coefficients, and Q{x,dx) and P{x,dx) are Z{H)- 
valued differential operators such that 

(4.7) L{x, dx) = Q{x, dx)P{x, dx), 

then 

4>(x, z) := P(x, dx)'I>H{x, z) 

is a generalized bispectral function, satisfying a spectral differential equation in x and a 
spectral differential-translation equation in z. 

More precisely, if 

L{x,dx) = Q'{x,dx)g{x)~^P'{x,dx) 

for Z{H)-valued differential operators with quasipolynomial coefficients Q'{x, dx), P'{x, dx) 
and a nondegenerate Z{H)-valued quasipolynomial function g{x), then the function 
^{x,z) := P'{x,dx)e^^In satisfies 

[P'{x,dx)Q'{x,dx)g{x)~^) ^{x,z) = ^{x,z){hL){z), 

g{x)^{x,z) = 4>(x,z) {{bL){z)~^h{Q')b{P')) . 

The first equation is a spectral differential equation in x. The second one is a spec¬ 
tral differential-translation equation in z; the point here is that h{Q') and h{P') are 
differential-translation operators. 

The classification of the factorizations of the form ()4.7p for monic differential operators 
L, Q and P is simpler than the classification in Theorem l4.3l The point is that the adelic 
type condition is not needed in this case since the operators P{x, dx) and Q{x, dx) are not 
required to have rational coefficients. (For each such factorization P{x, dx) and Q{x, dx) 
have quasirational coefficients, recall Lemma 02] (a).) Since L(x,dx) and Q{x,dx) are 
Z(77)-valued, their vector kernels are nondegenerate 77-invariant subspaces of QVn, 
recall (j4.3l) . One easily shows that every nondegenerate 77-invariant subspace of QPn 
of dimension divisible by n is the vector kernel of a unique monic differential operator 
P{x,dx) satisfying (14.7p for some differential operators L{x,dx) and Q{x,dx) with the 
stated properties. Thus we obtain the following: 

The generalized bispectral Darboux transformations in Theorem If.'A for monic differ¬ 
ential operators L, Q, and P are classified by the nondegenerate H-invariant subspaces 
of QVn of dimension divisible by n. For such a sub space, one constructs the unique 
monic differential operator P{x,dx) with that vector kernel and the corresponding gen¬ 
eralized bispectral function ^{x,z); this differential operator P{x,dx) is Z{H)-valued 
with quasipolynomial coefficients. 

5. Bispectral Darboux transformations from matrix Airy functions 

In this section we give an explicit classification of all bispectral Darboux transfor¬ 
mations from the matrix Airy functions. We often make use of the spectral theory of 
matrix polynomials reviewed in the appendix. 

5.1. We recall the setting of the (generalized) scalar Airy bispectral functions from [2]. 
Fix complex numbers ao 0, a 2 ,..., o:n-i- Consider the scalar differential operator 

N-l 

MAi(x,dx) ■■= dx + '^ aidx~" + oqx, 
i=l 
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called a (generalized) Airy operator. For N = 2 this is the classical Airy operator up 
to normalization. It is well known that the operator Mx\{x,dx) has A^-dimensional 
vector kernel over C consisting of entire functions. Choose a basis of the kernel, 
1/^1 (x),..., For each function '^(x) € vker M, define 

(5.1) z) := ijj{x + z). 

The latter satisfies the bispectral equations 

MAi(x,9a;)V’ = zip and M\i{z,dz)tp = xip. 

The bispectral Darboux transformations from it (in the sense of Theorem 12.II) were clas¬ 
sified in [2] and played an important role in the overall classification of scalar bispectral 
operators. 

Now we move to the matrix situation. As in the matrix rank 1 case from Section IU 
we take = Q 2 ■= C. Define the function 

^'Ai(a;,^;) = ip{x,z)In 

on C X C, where is the identity matrix. Let Bi and B2 be both equal to the algebra 
of matrix-valued differential operators with polynomial coefficients. We will use the 
variables x and z for the first and second algebra, respectively. Let b: Bi ^ B 2 he the 
isomorphism given by 

b{x):=MAi{z,dz), b{dx):=-dz, 6(IF) := IF, VW € M„(C). 

The first two equations imply that 6(MAi(x, dx)) = 2 ;. We are in the bispectral situation 
of Proposition 12.21 

P{x,dx)^Ai{x,z) = ^Aiix,z)(bP){z,dz), \/P{x,dx) G Bi. 

Let Ki and K 2 be the algebras of matrix-valued polynomial functions in x and z, 
respectively. Denote by Ai and A 2 the algebras consisting of operators of the form 

q{MAi{x,dx)) and q{MA\{z,dz)) 

for matrix-valued polynomials q{t). Then 

b{Ai) = K2, b{K2) = Ai. 

The regular elements of the algebras Ki and K 2 are the matrix-valued polynomial 
function 9(x) and f(z) that are nondegenerate in the sense that their determinants 
are nonzero polynomials. By standard commutation arguments for matrix differential 
operators the sets r{Ki) of regular elements are Ore subsets of B^. The localizations 
Bi[r{Ki)~^ and B 2 [r{K 2 )~^] are isomorphic to the algebras of differential operators 
with rational coefficients in the variables x and z. 

Theorem 12.11 implies the following: 

Theorem 5.1. Let q{t) be a matrix-valued polynomial. Let Q{x.,dx) and P{x,dx) be 
matrix-valued differential operators with rational coefficients, such that 

(5.2) q{MAi{x, dx)) = Qix, dx)P{x, dx). 

Then 

$(x,z) := P{x,dx)^Ai{x,z) 
is a matrix bispectral function. More precisely, if 

L{x,dx) = Q' {x,dx)g{x)~^ P' {x,dx) 


NONCOMMUTATIVE BISPECTRAL DARBOUX TRANSFORMATIONS 


19 


for matrix-valued differential operators with polynomial coefficients Q'{x,dx), P'{x,dx) 
and a nondegenerate matrix-valued polynomial g{x), then $(x,z) := P'{x, dx)^ a\{x■, z) 
satisfies 

{P'{x,dx)Q'{x,dx)g{x)~^) ^{x,z) = ^{x,z)q{z), 

g{x)^{x,z) = ^{x,z) {q{z)~'^{bQ'){z,d^){bP'){z,d^)) . 

In Theorem 15.41 below we classify all bispectral Darboux functions ^{x,z) that enter 
in Theorem 15.11 corresponding to q{t) and P{x,dx) with nondegenerate leading terms. 

5.2. First, we describe the kernels of the differential operators L{x, dx) in Theorem l5.ll 
In the notation (|5.ip . 

A) = {dlfji{x, z))\^=x = dif;i{x,X), VA € C, 1 < z < N,j G N. 

Define the space 

QAn:= 0 #(x,A)CT 

AgC,l<i<AJgN 

Denote the standard basis of C"" by {ei,..., e^}. We have 

{MAi{x,dx) - A) V'p^XjA) = V^p“^^(x,A), 

and 

{'ijj^^\x, X)\l<i<N,0<j<k — 1} 

forms a basis of the kernel of {MAi{x,dx) — A)^ for all A € C and k € Z_|_, see [2]. One 
easily obtains from this that 

(5.3) X)ei I A G C, 1 < i < A^, j G N, 1 < / < n} is a basis of QAn- 
The fact that 'fi{x) are in the kernel of MAi{x,dx) implies that 

QAn= 0 #(x,A)C[x]" 

AgC,l<i<Af,0<i<A 

and that 

(5.4) {x, X)x’^ei I A G C, 1 < z < iV, 0 < j < fV - 1, A: G N, 1 < / < n} 

is a basis of QAn- 

Fix a monic matrix-valued polynomial 

d 

= Z] QjP, where aj G Mn{C), ad = In- 
j=0 

Proposition 5.2. For each matrix polynomial q{t) as above, the vector kernel of the 
differential operator q{MAi{x,dx)) over C is ndN-dimensional and consists of entire 
functions. More precisely, for each root X of det{q{t)) and a choice of Jordan chains of 
q{t) corresponding to X as in Theorem \7. 1\ jvp.i, - - - Vk,.i | 1 < ^ < s}, the elements 

l<i<N,0<j <ki 

r=0 

belong to vker q(MAi(x, clx))- The set of all such elements for all roots X of det{q{t)) and 
a corresponding set of Jordan chains as in Theorem \ 7. l\ forms a basis of vker q{MAi{x, dx))- 
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Proof. Using the bispectrality of the functions 'ilji{x,z) and the definition of Jordan 
chains (IZ21), we obtain 


MAi{x,dx) A) 




^ 5^(MAi(x, z))\z=x 


Kr=0 


r=0 
j r 


t\ 


^^ii^pi{x, z)q{z))\,=x^^^ = = 

r=0 r=0 c=0 ^ ^ 


E 


c=0 




c: 


E 


(r - c)! 



To prove that the set in Proposition 15.21 gives a basis of vker ( 7 (MAi(x, (9a;)), first note 
that by ra the cardinality of this set equals Ndeg(det(q(t))) = ndN. Because 
dimvker g(MAi(x, 5a;)) E ndN, all we need to show is that these elements are linearly 
independent. This follows from (|5.3n and the linear independence in Theorem 17.11 □ 


A direct argument with matrix polynomials leads to the following corollary of Propo¬ 
sition 15.21 

Corollary 5.3. (a) Let q{t) be a monic matrix polynomial. If Q{x,dx) and P{x,dx) 
are monic matrix-valued differential operators with meromorphic coefficients in C which 
satisfy (15.2h . then vker P{x,dx) is a nondegenrate subspace of QAn. 

(b) For every nondenerate finite dimensional subspace V of QAn whose dimension 
is a multiple of nN, there exist a monic matrix polynomial q{t) and monic matrix¬ 
valued differential operators Q and P satisfying the properties in part (a) such that 
vker P(x, dx) = V. 

Introduce a ZA-action on QAn, where the generator cr of Za acts by 

a{ifil^\x, X)v) := X)v, V j € N, 1 < i < A7, u € C" 

and the i + 1 term in the right hand side is taken modulo N (i.e., N -\-l := 1). Here we 
use the fact that the set in (j5.3p is a basis of QAn. In terms of the basis elements from 
(15.4p . this action is given by 

a['ip\^\x, X)p{x)) := X)p{x), ^ 0 < j < N — 1,1 < i < N,p{x) € C[x]". 

This follows at once from the fact that the functions fifix) are in the kernel of M^fix, dx). 

The next theorem is our classification result for all bispectral Darboux transformations 
from the matrix Airy functions. 

Theorem 5.4. Assume that q{t) is a monic matrix-valued polynomial, and Q{x, dx) 
and P{x, dx) are monic matrix-valued differential operators whose coefficients are ratio¬ 
nal functions. If (j5.2p is satisfied, then vkei: P{x,dx) is a nondegenerate Zn- invariant 
subspace of QAn that has a basis consisting of elements of the form 

N-l 

(5.5) '^fi^^\x,X)pj{x), 

j=0 

where Pj{x) G C[x]"’. The indices i = 1,... ,N and the complex numbers A G C can be 
different for different basis elements. 
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Furthermore, for every such subspace of QAn of dimension that is a multiple of nN, 
there exists a matrix polynomial q{t) and a pair of operators Q and P with the above 
properties such that vker P{x,dx) equals this subspace. 

Corollary 5.5. The bispectral Darboux transformations from the matrix Airy function 
as in Theorem \5.1\ for a(t) and P{x,dx) with nondegenerate leading terms are 
exactly the functions of the form 

z) = W[x)P{x, dx)'^Ai{x, z), 

where P{x, dx) is a monic matrix differential operator of order dN with ndN-dimensional 
vector kernel which is a Tja- in-variant subspace of QAn possessing a basis consisting of 
elements of the form ([53]) and W (x) is a matrix-valued nondegenerate rational function. 
The operator P{x,dx) is uniquely reconstructed from the basis using the second part of 
Proposition \3.4\ 

5.3. We proceed with the proof of Theorem 15.41 The first part of the theorem follows 
from the following proposition and Corollary 15.31 ia). 

Proposition 5.6. Assume that P{x,dx) is a matrix-valued differential operator with 
rational coefficients. If 

N N-1 

p.6) eee X)pj fx) € vkei P{x,dx) 

AeC i=l j=o 

for some vector-valued polynomials pj^fx), then 
(n-i \ 

I 'f\^\x,X)p)i{x)\ ^ ykei P{x,dx), V A E C, 1 < * < 0 < /c < — 1. 

I- ’ / 

Proof. It follows from (j5.4p that 

A)e; |AeC, 0<j<iV — 

is a linearly independent set over C(x). Thus, (15.6p and the assumption that P{x,dx) 
has rational coefficients imply that 
N-l 

(5.7) X)pji{x) € vker P(x,dx), VAeC, l<z<Ai. 

j=0 

Using the facts that each such element belongs to vker P(x,dx) and P{x,dx) has rational 
coefficients, and recursively expressing 

N-l 

'tpi^\x, A) = - ^ X) -ao{x + X)'tfi{x, A), 

j=0 

produces a C[x]-linear combination of 

X)ei I 0 < j < — 1,1 < Z < n} 

that equals 0. In view of (|5.4p . this is only possible if each coefficient of the combination 
equals 0. However, exactly the same thing would happen if we apply P{x, dx) to 

N-l /N-1 

^)pii{x) = X] x)pYx) 

j=o y j=o 
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instead of (1521). Thus, for all 0 < /c < — 1, these elements belong to vker P{x,dx) 

which completes the proof of the proposition. □ 

The second part of Theorem 15.41 follows from the following proposition and Corollary 
E3](b). 

Proposition 5.7. Assume that V is a nondegenerate ndN-dimensional -invariant 
subspace of QAn having a basis consisting of elements of the form ()5.5I) . where Pj{x) € 
C[x]”. (The indices i = 1,... ,N and the complex numbers A € C can be different for 
different basis elements.) Let P{x,dx) be the unique monic matrix-valued differential 
operator on C of order dN whose vector kernel kernel equals V (recall Provosition \3.3\ 
(b)). Then P{x,dx) has rational coefficients. 

Proof. Denote the operator 

dN-l 

P{x,dx) = ai{x)dl, 

1=0 

where afx) are matrix-valued meromorphic functions on C. 

The definition of the Z^v-action on QAn implies that we can assume that the basis of 
V, consisting of elements of the form (15.5h . is itself Z^v-invariant. Let 

( N-l \ N-1 

Y = Y >')Pj{x), l<k<N-l 

j=0 ) j=0 

be N elements in the vector kernel of P{x,dx). Using the condition that they belong to 
vker P(x, cla;) and expressing 

Af-l 

= - Y 0ii'f^i\x,X) - aQ{x + \)fi{x,X), 

i=i 

produces a C[x]-linear combination of 

\)ei I 0 < j < — 1,1 < / < n} 

that equals 0. Once again ()5.4p implies that each coefficient of the combination should 
equal 0. Thus, the fact that one of the elements ()5.8p belongs to whei P{x,dx) is equiv¬ 
alent to imposing N conditions of the form 

dN-l 

Y^ ai{x)cij{x) = Q, V0<j<A^-l, 

1=0 

where cij{x) are some vector-valued polynomial functions uniquely determined from 
Pq{x), ... ,pN-iix). However, the condition that every element in ()5.8I) belongs to 
ykev P{x,dx) is equivalent to imposing exactly the same N conditions because iffxA) 
satisfy the same Airy equation. 

So, the conditions on the operator P{x,dx) (having vector kernel equal to V) are 
equivalent to imposing dN linear vector conditions on the matrix meromorphic functions 
{ai{x) I 0 < / < dN — 1}. Because the coefficients of these conditions are vector-valued 
polynomial functions and there is a unique monic differential operator P(x, dx) with the 
stated properties, all its coefficients afx) should be rational. □ 
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6. Examples 

In this section we present several examples illustrating the classification results from 
the previous two sections. 

6.1. First, we give three examples of matrix rank 1 bispectral Darboux transformations, 
classified in Theorem 14.31 and Corollary 14.41 

Example 6.1. In the case of 2 x 2 matrices, consider the differential operator 

L{x,dx) := dlh G Ai. 

(Below all instances of the identity matrix which are clear from the setting will be 
suppressed.) The vector kernel of L consists of all vector linear functions. In the setting 
of Theorem 14.31 take the subspace V of QV 2 with basis 

fi{x) := {x,0y, f 2 {x) = {a,xY 

for an arbitrary a G C. Set 

Fi{x) ■■= {fi{x),f 2 {x)) = (^ . 

Consider the corresponding operator P{x, dx) given by (13.31) . which in this case simplifies 
to 

P{x,dx) = dx-F[{x)F^\x) = dx-{^l^ 

By Theorem 14.3( we have the factorization L{x,dx) = Q{x,dx)P{x,dx) from which the 
operator Q is computed to be 

Q{x,dx) = dx+(^^~^' ““5 

This gives rise to the bispectral Darboux transformation 

^x,z) :=P{x,dx){e^^l2) = e^^ 

It satisfies 

L{x,dx)^{x,z) = ^{x,z)z‘^ 

where 

L{x, dx) ■■= Pix, dx)Qix, dx) = d‘^ - 

For the dual spectral equation, given by Theorem HU we need to represent L{x,dx) in 
the form (14.2h : 

d^h = Q'{x,dx)^P'{x,dx) with Q' := dxx‘^+(^ , P' := x‘^dx-(^ G Bi. 

Now, by Theorem 14.11 the function <h(x,z) satisfies the dual spectral equation 

<I>(x, z)K{z, dz) = z) 
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where 


K{z,dz) := 


dxX + 


1 


zd'^ -dz + 


X —a 
0 X 

0 —a 
0 0 


o b 


X^dr — 


X —a 
0 X 


diz + dz + 


0 a 
0 0 


= 2 \zdi-dz){di.z + d^)l2 + 


0 —2az ^02 
0 0 


When a = 1, this is the first example in [6] and [TH], where a (slightly more complicated) 
third order dual spectral equation was considered. 

Example 6.2. For 3x3 matrices, consider the operator L = G Ai. Its vector 
kernel consists of vector linear functions. Let V be the subspace of QV 3 with a basis 
consisting of the columns of the matrix 

/x 1 0^ 

Fi{x) := I 0 X 1 
yO 0 Xj 

The corresponding operator P{x,dx), given by (j3.3p . is 

P{x,dx) := dx - Fi{x)Fi{x)~^ = dx - {x~^l 3 - x~^J + 
where J is the super-diagonal matrix 



By Theorem 14.31 we have the factorization L(x,dx) = Q{x,dx)P{x,dx) from which the 
operator Q is computed to be 

Q{x, dx) = dx + {x~^h - x~‘^J + x~^j‘^). 

The bispectral Darboux transformation 


■.= P{x,dx){e^^h) = e^ 



— X 

x~‘^ 
z — X 


-3 


-1 


satisfies 

where 


L{x,dx)^{x, z) = ^{x,z)z‘^ 


L{x,dx) ■= P{x,dx)Q{x,dx) = d'^ — 2x '^Is + Ax ^J — 6x ^j‘^. 

To write the dual spectral equation in Theorem 14.11 we represent L{x,dx) in the form 
(02]), d^h = Q'{x,dx)x~^P'{x,dx) with 

Q'{x, dx) := dxX^ -b (x ^/3 — xJ + J^), P'{x, dx) '■= x^dx — {x^h — xJ + J^) G Bi. 

Theorem 14.11 implies that ^{x,z) satishes the dual spectral equation 

<h(x, z)A{z, dz) = x‘^^{x, z) 

with 

A{z,dz) := z~‘^b{Q'){z,dz)b{P'){z,dz) 

= z “^{zd^I^ — {dzis + dzJ + J‘^)){d^zl3 + {dz^s + dzJ + T^))- 
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0 

0 

1 

( 0 

0 

0 

0 

0 

0 92 + 


0 

0 

1 

0 

0/ ' 

[- 2 z -^ 

1 

0 


This bispectral function is the second example in [18], where a conjecture is stated about 
the full algebra of dual differential operators with eigenfunction d>(a:, z). It contains lower 
order operators, e.g., the degenerate operator 

1 0 0 \ 

0 0 I + I -2z-^ 2 0. 

0 0 1 / 

Example 6.3. For 3x3 matrices, consider the operator 

(dl dx 1 \ 

L{x,dx) = + dxJ + = I 0 I G Ai. 

\0 0 d^J 

A basis of vker L(x, dx) is given by 

(1, 0,0)^ (x, 0,0/, (0, 1, 0)^ (-xV2, X, 0)^ (0, -xV2, x)\ (-xV2, 0, 1)*. 

Consider the subspace of QV^ with a basis consisting of the columns of the matrix 
function 

/x —x^/2 

Fi(x) = j 0 X 0 

\0 0 1 ^ 

The corresponding operator P{x,dx), given by ()3.3I) . is P{x,dx) '■= dx — F{{x)Fi{x)~^, 
or explicitly, 

T-i Q/o ../o\ /i -2,x/2 -x^/2"' 


'x " —3/2 —x/2' 


P{x,dx) = dx- \ 0 X ^ 

0 0 


0 

0 


= X 


-{ 


xa,. - 10 1 0 

0 0 0 


} = x"^P'. 


By Theorem 14.31 we have the factorization L{x,dx) = Q{x,dx)P{x,dx) from which the 
operator Q is computed to be 

T-i I/O ^/o\ /i _xl2 -xV2^ 

= ^dxX + 0 


'X " —1/2 —x/2' 


Q = dx+ \ 0 X 
0 0 


-1 


1 

,0 0 


X 

0 


|x ^ = Q'x 


The operators P'{x, dx) and Q'{x, dx) will be needed for the dual bispectral equation. 
The bispectral Darboux transformation 

"-i 3/2 x/2^ 


$(x,z) ■.= P{x,dx){e^^l3) = 


z — X 

0 
0 


z — X 

0 


-1 


0 


satisfies 


_ / z^ z 1 

L{x,dx)^{x,z) = ^{x,z) \ Q z^ z 

\0 0 z^ 

where the differential operator L{x,dx) is given by 

L(x, dx) := P{x, dx)Q{x, dx), 

omitting the detailed expansion for the sake of brevity. By Theorem sn the function 
^>(x, z) satisfies the dual spectral equation 
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where 

A(z, dz) = {z^h + zJ + j‘^)~^h{Q'){z, d^)b{P'){z, 5^) = 


0 

0 


— 2 ; 

0 


-3 


-3 


— Z 

,v-2 


A d,/2 -dl/2^ 

- ( 0 1 -dz 

,0 0 0 


||5^z + 



6.2. We finish with an example illustrating the classification results for all bispectral 
Darboux transformations from the matrix Airy functions. 

Example 6.4. Consider the matrix polynomial 

q{t) = {ht - Jf = ht^ - 2Jt 

where 

J := 


For the second order Airy operator Mx\{x,dx) = d'^ — x, consider the operator 
L{x,dx) := q{MAi{x,d^)) = {hdl - hx - jf € Bi. 

It satisfies h{L{x,dx)) = q{z) = z^l 2 — 2zJ, i.e., 

L{x,dx)'i>Ai{x,z) = '^Aiix,z){z‘^l 2 - 2zJ). 

According to our scheme (from Proposition 15.21 Theorem 15.41 and Corollary 15.5p . we 
first need to find the Jordan chains of the monic matrix polynomial q{t). Its character¬ 
istic polynomial x(t) = det{q{t)) = has only one root 0 and 

^(0) = 0, g'(0) = -2J, q"{0) = 2 / 2 . 

Denote by 61,62 the standard basis of C^. Then 

62 , 

61 , ( 1 / 2 ) 62,61 

is a maximal family of Jordan chains as in Theorem l7.ll Recall from llS.ll that 'ijji{x),ip 2 {x) 
denotes a basis of vker MAi(x, cla,). Proposition 15.21 produces the following basis of 
vker L(x, (9a;) : 


i/i(x) + V’f(a^)/2 
V’'(x)/2 


i = 1,2. 


0 A (Mx)\ ( ^i{x) 

Mx)) ’ V 0 y ’ \Mx)/‘^, 

(One can further simplify it by taking linear combinations.) Consider the 4-dimensional 
subspace V of QA 2 with basis 


■(/'(x) J- ijji{x] 

ipiix) 


i/f(x) 

V('(x) 


i = 1,2. 


It satisfies the conditions in Theorem 15.41 in particular V is Z 2 -invariant. It can be 
computed by using quasideterminants. However, more effectively, it can be computed 
by setting P{x,dx) = 8“^ + h{x)dx + c(x) and then solving for h{x) and c(x) from the 
equations P{x,dx)f{x) = 0 for the 4 basis elements /(x) of V. The concrete solution 
of the equations for b{x) and c(x) is done by expressing all derivatives (x) in terms 
of tpi^x) and V’((x) using the Airy equation, and then setting the coefficients in front of 
'il>i{x) and i/((x) in the equations P{x,dx)f{x) = 0 to be equal to 0. This gives that 


P{x, dx) = dl + 


1 


X — 1 


0 1 -X 
1 -1 


dx - 


1 


X — 1 


(X-1)2 

0 


2x-2 
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By Theorem 14.31 L{x,dx) = Q{x,dx)P{x,dx) for some second order monic differential 
operator Q{x,dx)- By comparing coefficients one obtains 


Q(x, dx) 


di 


1 /O 1 — 1 /— 

VI -1 ) ^ V 0 




The bispectral Darboux transformation 


$(x,z) ■.= P{x,dx)'^A\{x.,z) 


satisfies 

L(x, dx)^{x, z) = ^>(x, z){z^l 2 — 2J) 

where 

L{x,dx) := P{x,dx)Q{x,dx). 

For the sake of brevity we leave the detailed computation of L{x,dx) to the reader. 

To write the dual spectral equation from Theorem 15.11 we represent q{Mf^\{x, dx)) in 
the form (|4.2|) . 

q{Mxi{x,dx)) = Q'{x,dx)-. — \^P'ix,dx) 

(X — 1 )^ 

with 

Q'{x,dx) := Q{x,dx){x - 1 ), P'{x,dx) := (x - l)P{x,dx) G Bi- 


Explicitly we have 


Ph t l-x\„ /(x - 1)^ 2x - 2\ 

P (x, = (x - 0 x2 )’ 

Q'{x,dx) = d^{x-l)-dx(j + ( 0 -(x-l)^)- 

Theorem 15.11 implies that ^{x,z) satisfies the dual spectral equation 


<h(x, z)k{z, dz) = (x — l)^<l>(x, z) 


with 


K{z,dz) := {zh - J)-%{Q'){z,dz)h{P'){z,dz) 

= {z~‘^ + 2z~^J)b{Q'){z,dz)b{P'){z,dz) 


where 


HP'){z,dz) 

b{Q')iz,dz) 


{Mz - 1)^2 - 

d^.iMz -1)^ + dz 


i-_Mzy^_^iMz-i)^ 

•0 T-Mi _ 


2Mz - 2\ 

Ml ) 

{M, - 1 )V 


and we abbreviated Mx := Mx\{x,dx)-, Mz := Mp^\{z,dz)- 
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7. Appendix. Matrix polynomials 

Here we collect some facts about the spectral theory of matrix polynomials, used in 
the paper. We follow |15] and refer the reader to m Sections 1.4-1.6 and SI.5] for 
details. 

Consider a monic matrix polynomial 


d 

ajt^, where Oj E Mn{C), ad = In- 

j=0 


Define the monic polynomial 

X(t) := det(q(t)) 

of order nd. It equals the characteristic polynomial of the companion matrix 


(7.1) 


C = 


( ° 

/ 

0 

0 \ 


0 

0 

I 

0 






€ Mnd{II)i 

\-ao 

-ai 

-02 

... —Ud-l) 


U0,Ui . 


E C”, 

with uo 7 ^ 0 

is called a Jordan chain of q{t) 


see m Theorem 1.1]. 

A sequence of vectc 
of length k + 1 corresponding to A E C if the following equalities holds: 


(7.2) = 0, VO < j < A:. 

r=0 

Such a chain exists only if A is a root of x(t). Its leading term vq is an eigenvector of 
q{X). The condition (17.2h is equivalent to saying that the function 

j,k—l \ 

is a solution of the ODE q{dx)u{x) = 0. When q{t) is the characteristic polynomial of a 
square matrix A, this notion recovers the classical notion of a Jordan chain for A. 

The following theorem is a generalization of the well known Jordan normal form 
theorem for square matrices, see |15l Theorem 1.12] for details. 

Theorem 7.1. Let q{t) be a monic matrix polynomial and \ be a root of x{t) of multi¬ 
plicity m. Then there exist positive integers ki,... ,ks and Jordan chains 

vo,i, ■ ■ € C 

(for 1 < I < s) of q{t) with respect to A such that {uo,i,..., uo,^} is a basis of ker q{X) 
and 


(7.3) {k\ -|- 1) {kg + 1) — m. 

The lengths of the Jordan chains in the theorem can be read off from the local Smith 
form of q{t) in terms of the partial multiplicities of q{t) at A, see [TSl Section SI.5]. 
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